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ABSTRACT 

This  thesis  presents  an  analysis  of  a  direct  current 
conduction  magnetohydrodynamic  generator.   Laminar,  in- 
compressible viscous  flow  is  considered.   Experimental 
evidence  is  offered  to  show  that  the  laminar  assumption 
is  valid  in  the  presence  of  strong  magnetic  fields. 


Efficiencies  are  calculated  for  generators  of  various 
length  and  aspect  ratio,  efficiency  curves  are  presented 
and  a  method  of  driving  the  flow  is  considered. 
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Chapter  I 

INTRODUCTION 
1.1  Background 

The  minimization  of  sound  transmitted  from  a  submarine 
through  the  sea  to  the  listening  devices  of  a  potential 
enemy  is  a  major  problem  in  submarine  design.   A  large 
fraction  of  the  radiated  noise  is  generated  by  the  submar- 
ine's propulsion  plant:   the  main  engines  and  the  screws 
and  shafting 

Consequently  it  would  by  desirable  to  devise  a  propul- 
sion system  having  neither  shaft  nor  screw  and  a  minimum 
quantity  of  rotating  machinery.   Doragh*  (l)  has  shown  the 
feasibility  of  propelling  a  submarine  by  means  of  a  sea 
water  jet  accelerated  in  a  magnetohydrodynamic  channel. 
This  thesis  will  be  a  study  of  the  closely  related  problem 
of  providing  the  electrical  energy  necessary  to  drive  the 
magnetohydrodynamic  flow. 

The  pioneer  figure  in  the  study  of  magnetohydrodynamic 
channel  flow  was  Hartmann  (2).   In  1937  he  made  the  first 
detailed  theoretical  analysis  of  the  flow  of  a  conducting 
fluid  in  a  high  aspect  ratio  channel  under  the  influence 
of  a  uniform  transverse  magnetic  field.   Hartmann  and 
Lazarus  (3),  in  that  same  year,  reported  on  observations 
made  on  mercury  flow  in  a  transverse  magnetic  field.   To 
date  the  principal  use  of  magnetohydrodynamic  channel  flow 
has  been  as  electromagnetic  pumps  to  circulate  liquid 
metal  coolants  in  atomic  power  installations. 

At  the  present  time  there  is  considerable  interest  in 
the  development  of  magnetohydrodynamic  power  generators. 


*  The  bracketed  numbers  refer  to  the  Bibliography  which 
begins  on  Page  28 


2. 

Most  of  these  proposed  generators  make  use  of  a  conducting 
gas  or  plasma  as  the  working  fluid.   This  type  of  magneto- 
hydrodynamic  generating  system  shows  promise  of  supplying 
large  quantities  of  electric  power  at  thermal  efficiencies 
unattainable  with  present  equipment  (L|.,5,6,7)  • 

1. 2  The  Object  and  Scope  of  This  Study 

The  object  of  this  study  is  the  design  of  a  liquid 
metal  magnetohydrodynamic  direct  current  conduction  gener- 
ator compatible  with  and  utilizing  the  same  magnetic  field 
as  the  magnetohydrodynamic  pump  designed  by  Doragh.   As 
will  be  shown,  the  use  of  a  liquid  metal,  instead  of  a  con- 
ducting gas  or  plasma,  as  the  working  fluid  permits  the 
design  of  a  complete  energy  conversion  cycle  devoid  of  mech- 
anical pumps  and  compressors.   Since  the  basic  goal  of  this 
design  is  to  achieve  a  system  free  of  these  mechanical  noise 
sources,  gaseous  working  fluids  will  not  be  considered. 

The  channel  under  consideration  is  rectangular  in  cross 
section  and  of  high  aspect  ratio.   The  applied  transverse 
magnetic  field  is  uniform  and  parallel  to  the  shorter  sides 
of  the  channel  cross  section.   Furthermore,   the  liquid 
metal  is  assumed  to  be  isotropic,  homogeneous  and  incom- 
pressible 

The  study  will  be  restricted  to  laminar  viscous  flow, 
since  there  is  little  or  no  experimental  data  available 
concerning  turbulent  flow  of  an  electrically  conducting 
fluid  in  a  strong  transverse  magnetic  field-   In  fact,  the 
assumption  of  laminar  flow  may  be  quite  valid  since  there 
is  some  experimental  evidence  that  indicates  a  transverse 
magnetic  field  delays  the  transition  from  laminar  to  turb- 
ulent flow.   Murgatroyd  (8)  in  his  experiments  with  mercury 
found  that  laminar  flow  was  stable  when  the  hydraulic 
Reynolds  number  was  less  than  225  times  the  Hartmann  number. 
Lock  (9)  hypothesized  that  laminar  flow  would  be  stable  for 
a  Reynolds  number  less  than  lj.2,000  times  the  Hartmann  number 


3. 
for  larger  Hartmann  numbers.  Since  this  study  will  Toe  deal- 
ing only  with  large  Hartmann  numbers,  both  theory  and  exper- 
iment predict  laminar  flow. 
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Chapter  II 
PROCEDURE 


2.1  Basic  Equations 


The  fluid  under  consideration  is  homogeneous,  isotropic, 
non-magnetizable,  non-polarizable  and  incompressible   It  is 
characterized  by  a  permeability  p,0,  an  electrical  conductivity 
a  mass  density  p  and  a  dynamic  viscosity  T\  .      The  equations 
that  govern  the  motion  of  such  a  fluid  are  the  magnetohydro- 
dynamic  Navier-Stokes  equation. 

p^Z.  =,  -grad  p  +  v^v  +  JXB,  (2.1.1) 

J./ T/ 

the  incompressibility  condition 

div  v  =  o,  (2.1. 2) 

Maxwell's  equations 

curl  H  =  J,  (2  1.3) 

curl  E  =  -|*^f,  (2.14) 

and 

div  H  =  o;  (2.1.5) 

and  the  constitutive  relation 

J  =  ct(e  +  v  X  |iH).  (2.1.6) 

Where  H  is  the  magnetic  field  intensity,  J  is  the  conduction 
current  density,  E  is  the  electric  field  intensity,  v  is  the 
fluid  velocity  and  P  is  the  pressure.   Rationalized  MKS  units 
are  used  in  all  equations  and  all  equations  are  written  with 
respect  to  a  fixed  coordinate  system. 


2. 2  Pressure  and  Velocity  Relations 

The  geometry  of  the  generator  to  be  considered  is  il- 
lustrated in  Figure  I,   Fluid  flows  in  the  z -direction 
with  a  velocity  that  is  a  function  of  both  x  and  y.   A 
uniform  magnetic  field  H  is  applied  in  the  x  direction. 
The  channel  walls  at  y  =  -  w/2  are  highly  conducting 
buses  that  are  maintained  at  a  potential  difference  by 
the  external  constraint  imposed  by  the  load  resistance. 
The  walls  at  x  =  -  d/2  are  non-conducting.   The  electric 
field  E  resulting  from  the  transverse  potential  gradient 
is  assumed  to  be  in  the  negative  y  direction,  and  the 
current  density  J  is  assumed  to  be  in  the  y  direction. 

The  channel  cross  section  is  shown  in  Figure  II. 
Load  current  I  is  assumed  to  flow  symmetrically  in  the 
external  circuit   This  insures  that  the  magnetic  field 
resulting  from  the  flow  of  load  current  will  be  in  the 
z  direction  and  so  will  not  affect  the  generator  oper- 
ation.  End  effects  will  not  be  considered  at  this  time. 

Under  these  assumptions,  the  z  component  of  Equation 
(2.1.1)  becomes 

0  =  -*P  +^v2   +  ^v-   -J  B  (2XD 


&X*      '& 


The  elimination  of  J  by  use  of  Equation  (2.1.6)  yields 


*2 


V 


.2 


+  Vv   -  **B2v  +  <T  (EB  -_L  ^  P)  =  0  *     (2.2.2) 


bx    *y2   ^      1)      ff  * 

This  equation  can  be  written  as 


*  From  this  point  on  indices  indicating  vector  directions  will 
be  omitted.   All  directions  are  as  defined  in  Figure  I. 


6. 


B 


f  A  V 


X 


z 


y 


j 


Figure  I 
Channel   Geometry 
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Figure  II 
Channel  Cross  Section  and  External  Electric  Circuit 


> v     a 


+ 


'v      ,2 


^x2   £  y 


^-  -  k  v  +  P  =  0  (2.2.3) 


where 

k2  =  B2  (d/2)2     =  M2 

-(d/2)2  (d/2)2 


(2.2.10 


and 


M  =  B  (d/2)  ^21  (Hartmann  number)  (2.2.5) 


P  =  51(EB  -  i^JP)  (2.2.6) 

*1      ^  *  z 


A  solution  to  Equation  (2.2.3)  of  the  form 

2(x)    cos    (2n   -   pity        dy  {22  ?) 


2n 


is   desired,   where 

v2n-l<x)   =  -     \     v  cos    (2n"l)lty       dy  (2.2.8) 


W      '-XL  w 


If  Equation  (2.2.3)  is  multiplied  by  the  function 
-  cos  \dn~    )J     ^y  ancj  the  resulting  equation  integrated 


w 

w 


over  the  interval  y  -  -w/2  to  y  =  w/2,  and  the  boundary 
condition,  v(x,  -w/2)  =  0  is  applied,  the  result  is 


~1 5 \     f*/o\2  +     L        w         J  Un-l)n: 


b   ^ 


or 


8. 
,n+l 

(2.2.9) 
>Sa-l   "  Mn  v2n-l   +  ^"'^  "  °  (2-2.10) 


>x2  (2n-l)n 

where 

M  2     =  M2  +    t2""1)* 

n        13727  2     —  <2-2-") 


Define  the  function  , 

°2n-l     n  v2n-l   U-2.J^; 

(2n-l)7c 

then  Equation  (2.2.10)  can  be  expressed  as 


^2n-l  -  Mn2£  ^^  -  0  (2.2.13) 

^x2 


This  has  the  solution 

62n_1  =  Acosh  Mnx  +  B  sinh  Mnx  (2.2.34) 

where  A  and  B  are  arbitrary  constants  to  be  evaluated  from 
the  boundary  conditions. 

Because  of  the  symmetry  of  the  problem,  v,  and  hence 
6p  -.  must  be  an  even  function  of  x.   Therefore  B  =  0. 
A  may  be  evaluated  from  the  requirement  that  v(-d/2,y)  =  0. 
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The  solution  for  6„  .  is  then 

2n-l 


2n-l 


-^P(-l) 


n  +  1 


M  d 


(2n-l)rc  cosh  n 


(2.2,15) 


and 


_  l4-P(-l)n  +  -1   -,    cosh  Mx 
v2n-l   -1   2_ 


(2n-l)nM 


n 


whence 


c^ 


v  =  ^P 


% 


(-D 


n  +  1 


v^~\ 


(2n-l)Mn' 


cosh  M  d 
n 


(2  2.16) 


t   cosh  M  x 
1-       n 


cosh  M  d 


cos 


(2n-l)l  n  y 


w 


(2.2.17) 


The  expression  for  the  mean  velocity  y  can  be  found  by 
integrating  Equation  (2.2.17)  across  the  channel  with  the 
result 

Crt 


V    = 


8P 


k2     £rr      (2n-l)2Mn2 


1    -  ±_  tanh    (Mn|) 


Mnd 


(2.2.18) 


Equations  (2.2.6)  and  (2.2.18)  can  now  be  combined  to 
give  the  pressure  gradient  down  the  channel. 


10, 


Po  =  *  P   =(TEB  -  ^ 


(2n-l)T! 


1  -  f tanh(M  f) 

2„2L   *?      ^ 


o^  N    '  n 


(2.2.19) 


2.3  Terminal  Impedance 

For  large  Hartmann  numbers  Equation  co2.2.17)  predicts 
that  the  normal  parabolic  velicoty  profile  associated  with 
laminar  flow  is  distorted  by  the  presence  of  the  transverse 
magnetic  field.   The  flow  in  the  center  of  the  channel  is 
retarded  and  that  near  the  walls  accelerated  so  that  the 
result  is  a  core  of  uniform  flow  and  a  boundary  layer 
whose  thickness  is  of  the  ordere  of  d/M.  ** 

In  determining  the  terminal  characteristics  of  the 
generator  the  velocity  will  be  assumed  to  be  uniform  and 
equal  to  v_   The  Hartmann  number  will  be  large  enough  to 
make  this  assumption  valid. 

With  the  assumption  of  a  uniform  velocity  profile, 
the  equivalent  electrical  circuit  for  the  generator  can 
be  determined  quite  easily  from  Equation  (2.1.6).   After 
transforming  this  equation  from  field  variables  to  cir- 
cuit variables  the  result  is 

lid  .  (vB  +  I).  (2.3.1) 

w 

I  and  V  are  positive  in  the  directions  defined  in  Figure  II 

*  For  a  complete  discussion  of  these  results  see  Ryabinin 

and  Khozhainov  (10). 

**  See  Shercliff  (11 )  page  137- 
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This  current-voltage  relationship  determines  the  equivalent 
electrical  circuit  shown  schematically  in  Figure  III. 
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Figure   III 

A  characteristic  of  this  generator  is  its  low  terminal 
voltage.   This  is  readily  seen  from  the  expression  for  open 
circuit  voltage,  vBw.   For  a  liquid  metal  the  maximum  mean 
velocity  practiable  is  of  the  order  of  10  meters  per  second. 
The  channel  width  is  of  the  order  of  one  meter.   The  B  field 
used  by  Doragh,  and  hence  the  field  acting  on  this  generator 
is  ten  Webers  per  square  meter.   The  open  circuit  voltage  then 
is  limited  to  a  maximum  of  about  100  volts. 
2.1+  Power  Relations 

Define  R. ,  the  internal  resistance  of  the  generator,  by 


Ri  = 


w 


<*ld 


The  terminal  voltage,  V,  can  then  be  expressed  as 


(2.^.1) 


,,   vBwRT 

v  .  ~ Jj , 

Ri  +  RL 


(2.I+.2) 


and  the  load  current,  I,  as 
T  _  vBw 


Ri  +  RL 


(2-it-3) 


12, 


The  electric  power  delivered  to  the  load,  P,  ,  is  given  by 

(R±   +  RL)2 


And  the  generator  efficiency,  e,  is  given  by 

e=^ 

Rate  of  doing  mechanical  work  on  the  fluid 

P  v  lwd       (R.  +  RT)2P  v  lwd 

O—  x  1      L'   O — 

2.5  The  Energy  Conversion  Cycle 

The  source  of  energy  for  operating  the  generator  is 
assumed  to  be  a  nuclear  reactor.  The  problem  of  con- 
verting the  thermal  energy  of  thw  reactor  into  mechani- 
cal energy  necessary  to  drive  the  generator  without  the 
intervention  of  mechanical  pumps,  compressors  and  other 
mechanical  noise  sources  is  a  difficult  one.  The  cycle 
considered  here  was  devised  by  Elliot  (12).  It  is  shown 
schematically  in  Figure  IV 

Elliot  carries  out  the  thermal  to  mechanical  energy 
conversion  by  mixing  a  hot  liquid  metal  of  low  vapor 
pressure  with  a  cooler  liquid  metal  of  higher  vapor 
pressure.   The  temperature  of  the  hot  metal  is  above 
the  boiling  point  of  the  cooler  metal,  and  in  the  mixing 
process  the  lower  temperature  metal  vaporizes.   The  mix- 
ture of  vapor  and  liquid  metal  droplets  is  passed  into  a 
nozzle  where  the  vapor  expands  and  accelerates  the  liquid 
phase.   The  vapor  and  liquid  phases  are  then  separated 
from  each  other.   The  vapor  is  condensed  and  pumped  back 
to  tht  mixing  chamber  by  means  of  a  magnetohydrodynamic 


13 
pump.   The  liquid  flows  through  the  generator,  to  the  heat 
source  and  back  to  the  mixer. 

The  essential  feature  of  this  cycle  is  that  it  has  no 
moving  parts  and  should  satisfy  the  requirement  for  noise- 
free  operation  to  a  high  degree.   Its  major  drawbacks  are 
that  it  is  relatively  inefficient  and  the  problem  of  com- 
plete separation  of  the  two  phases  has  not  yet  been 
solved. 
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Figure  IV  A  Proposed  Power  Generation  Cycle 
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Chapter  III 

RESULTS 

3  1  Results 

Generator  efficiencies  were  calculated  for  generators 
of  several  lengths  and  aspect  ratios  for  different  mean 
flow  velocities.   The  resultant  plots  of  efficiency  versus 
velocity  are  shown  in  Figures  V,  VI  and  VII.   Figures  VII, 
IX  and  X  show  characteristics  of  a  generator  one  meter 
wide,  one  meter  long  and  with  an  aspect  ratio  of  five. 

Liquidlithium  was  assumed  to  be  the  working  fluid  in 
all  cases. 
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Figure  X 
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Chapter  IV 
DISCUSSION  OF  RESULTS 

Lj..  1  Discussion 

The  results  show  that,  with  the  exception  of  the  thin 
high  aspect  ratio  channel  there  was  relatively  little 
change  in  efficiency  as  the  velocity  was  varied  through- 
out the  useful  range.,   Also,  the  longer  channels  were 
more  efficient  than  the  short  channels   This  result 
was  to  be  expected  since  the  very  low  load  impedence 
coupled  to  the  generator  caused  large  current  flows. 
These  currents  were  large  enough  that  the  viscous  los- 
ses in  the  fluid  were  insignificant  when  compared  with 
ohmic  losses.   Viscous  forces  were  important  only  in 
the  thin  channel  at  high  speeds. 

The  results  also  demonstrate  the  obvious  fact  that 
channel  efficiency  increases  as  channel  cross  section  is 
increased   Viscous  losses  are  basically  surface  effects 
and  anything  that  is  done  to  decrease  the  ratio  of  sur- 
face area  to  volume  of  the  channel  serves  to  lower  the 
viscous  losses. 

Ohmic  losses,  on  the  other  hand,  are  increased  as 
channel  width  is  increased  and  are  lowered  when  length 
or  height  is  increased 

At  this  time  it  might  be  worthwhile  to  test  the 
validity  of  the  assumption  of  laminar  flow  upon  which 
the  whole  analysis  was  based  For  liquid  lithium*  in  a 
magnetic  fluz  density  of  ten  webers  per  square  meter 
the  Hartmann  numbers  ranged  from  6  x  1J4"  to  2  x  10*^ 
The  hydraulic  Reynolds  number  for  these  same  flows 
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ranged  from  10J  to  10'    The  ratio  of  Reynolds  number 
to  Hartmann  number  was  at  all  times  less  than  the 
critical  value  of  225  found  by  Murgatroyd.   There  is 


*  For  the  properties  of  lithium  see  Appendix  A 
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still  some  question  as  to  the  validity  of  the  relationship 

for  the  critical  Reynolds  number  being  equal  to  225  times 

the  Hartmann  number  for  the  extremely  high  Hartmann  numbers 

encountered  here.   In  this  respect  it  is  interesting  to 

note  that,  in  his  experiments,  Murgatroyd  was  unable  to  ob- 

5 
tain  turbulent  flow  for  a  Reynolds  number  of  10  .. 

One  of  the  major  disadvantages  of  the  constant  velocity 
channel  is  the  large  pressure  gradient  required  to  produce 
any  significant  amount  of  electric  energy.   Since  the  liquid 
metal  is  incompressible  there  is  no  conversion  of  thermal 
energy  to  electrical  energy  within  the  channel .   The  kin- 
etic energy  of  the  flow  also  is  constant  in  the  channel 
Thus  the  difference  between  generator  inlet  pressure  and 
outlet  pressure  alone  must  account  for  both  channel  losses 
and  power  delivered  to  the  load. 

Another  disadvantage  of  the  liquid  metal  channel  as  a 
generator  is  the  low  terminal  voltage  available.   It  is 
desirable  to  operate  at  low  flow  velocities  to  keep  viscous 
losses  low.   Since  the  terminal  voltage  of  the  device  is 
dictated  by  the  product  of  magnetic  flux  density,  velocity 
and  channel  width,  a  terminal  voltage  of  no  more  than  one 
hundred  volts  is  expected,   This  low  voltage  requires  a 
low  load  impedence  in  order  to  extract  any  significant 
amount  of  electrical  energy  from  the  flow. 
L)..  2  Conclusion 

The  liquid  metal  conduction  generator  is  extremely 
attractive  because  of  its  small  size  and  simplicity.   The 
major  disadvantages  listed  above  are  serious  but  not  in 
themselves  sufficient  grounds  for  discarding  the  device. 

It  is  only  when  one  looks  further  than  the  generator 
itself  that  insurmountable  difficulties  appear.   The 
absence  of  an  efficient  method  for  converting  thermal 
energy  to  flow  energy  without  using  mechanical  equip- 
ment is  a  severe  handicap.   Elliot's  cycle  mentioned 
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earlier  has  many  drawbacks.   It  is  inherently  inefficient 
and  its  temperature  requirements  exceed  the  capabilities 
of  presently  available  materials.   The  cycle  requires  the 
use  of  a  liquid  metal  with  a  high  boiling  point  in  the 
liquid  loop.   It  follows  that  any  metal  meeting  this  re- 
quirement would  be  solid  at  room  temperature.   As  a  result 
a  means  must  be  provided  for  keeping  the  metal  from  freez- 
ing when  the  system  is  not  in  use.   If  lithium  is  used 
this  means  maintaining  a  temperature  of  above  35^4-°  F- 

The  flow  could  be  driven  by  a  mechanical  pump  opera- 
ting in  a  conventional  steam  cycle.   However,  this  would 
defeat  the  major  reason  for  having  an  MHD  generator: 
noise  elimination.   Besides,  if  a  conventional  steam 
cycle  must  be  built  to  operate  a  pump,  it  may  just  as  well 
drive  a  rotating  generator  directly. 
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Chapter  V 

CONCLUSIONS  AND  RECOMMENDATIONS 

5- 1  Advantages 

The  prime  advantage  of  the  generator  studied  is  its  small 
size  and  simplicity.   In  addition,  it  meets  the  basic  require- 
ment of  relatively  noise-free  operation. 
5. 2  Disadvantages 

The  disadvantages  of  this  type  of  generator  system  are 
manifold.   Terminal  voltage  is  low,  losses  in  the  system 
excluding  the  generator  are  relatively  high  and  there  is  a 
serious  materials  problem. 
5. 3 ■ Recommendations 

Since  the  difficulties  of  producing  electrical  energy 
by  means  of  a  liquid  metal  cycle  far  outweigh  the  advantages 
of  such  a  system,  as  study  of  an  MHD  power  generator  system 
using  an  ionized  gas  as  a  working  fluid  would  certainly  be  in 
order.   If  such  a  system  could  be  made  small  enough  in  phy- 
sical size  to  be  placed  in  a  submarine  it  might  very  well 
provide  the  quiet  operation  we  are  seeking. 


APPENDIX  A 
THE  PROPERTIES  OF  LIQUID  LITHIUM  (13) 

Melting  Point  35^°  F 

Boiling  Point  21+03°  F. 

Density 


Viscosity- 


Resistivity 
Hartmann  Number 


.507   gm/cc 

200°   C. 

.1+90 

1+00 

.1+71+ 

630 

.1+57 

800 

.iiia 

1000 

.5918 

Centerpoises 

183. 9° 

.571+9 

193.2 

•  551+1 

208.1 

1+917 

250.8 

.1+5^8 

205.5 

14-5.25 

micro  ohm-cm 

2300  c 

=  B(d)Jlr    = 
2    ^ 

6.82  x   105(d) 

C 
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APPENDIX  B 
METHODS  OF  COMPUTATION 

All  calculations  were  made  using  the  properties  of 
lithium  listed  in  Appendix  A 

Equation  (2.2.19)  was  solved  for  Po,  Equation  (2.2.I4.) 
for  PT  and  Equation  (2.I4..5)  then  solved  to  determine 
efficiency. 

A  magnetic  flux  density  of  ten  webers  per  square  meter 
was  assumed  in  all  calculations. 
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